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APPENDIX A 
EIGENSYSTEM SENSITIVITY MINIMIZATION AND CALCULATION 


A.1 EIGENVALUE DERIVATIVES 


Consider an nxn matrix A with n_ distinct eigenvalues 


Ay» A +, AW Let the diagonal matrix A_ be given as 


2? 


A = (A.1) 


Let the ith column of the matrix U,u; 


corresponding to the eigenvalue rGs Then 


be the right eigenvector 


AU = UA (%.2) 


Let the ith column of the matrix V,v; be the left eigenvector 
corresponding to the eigenvalue 5. Then 


H (A.3) 


< 
> 
u 
z 


or 


H H 


AV (A.4) 


u 
5 
u 
5 


where the "H'' superscript indicates the Hermitian transpose, 
or the complex conjugate of the ordinary transpose. 








— i 


When A _ is not symmetric, the right eigenvectors u; 


are not mutually orthogonal. However, the right and left eigen- 
systems have a property called biorthogonality. A right eigen- 
vector uj; and a left eigenvector vy corresponding to 


distinct eigenvalues AG and ae are orthogonal. 


Proof [52] 


By definition 


Au; = Aju; : Aveyis = Nav (A. 5) 





vjAu; = uae ; u;A we AiUsV, (A.6) 
The second relation gives 
vj Au; = AGV 5X5 (A.7) 


Subtracting this from the first inner product gives 
é = g 
(A; Aj) V5U; 0 (A. 8) 


If r; # Ais then Vj and u; are orthogonal. 


We can now show that the inverse of the matrix U isa 
matrix of left eigenvectors of A: 


cu-2)Ha = acuty# (A.9) 
Proof [52] 
Define V such that v" = p“!, By definition of U 
AU = UA (A.10) 
Z 





Since vily = uv = I, we have 


viau = A (A.11) 


via = avtl (A.12) 


Equation (A.12) is the same as Eq. (A.3). 


Now we Can consider eigenvalue and eigenvector derivatives. 
Differentiating Eq. (A.2) by some parameter © gives 


AgU + AUg = UgA + UAg (A.13) 
; -1 H 
Multiply through by U =V 
H H _ =H 

VAgU + AV Ug = V UgA + Mg (A.14) 


The i,j elements of Eq. (A.14) can be written as 


H H = H 
[V"AQU],; * [V'Ug];; a4 7 (V Uglis 44 * [Mg]jj (4-15) 
For i=j we have 
H vC 
[V AgU];; = 56 r; (A.16) 
or 
fo a0), , * ee (A.17) 
O° -ii 08 ‘i : 


A.1.1 Computation Using Real Arithmetic 


If the matrix A is not symmetric, then its eigenvalues 
and eigenvectors will be complex. During numerical computations, 
the use of complex arithmetic should be avoided; it is time- 
consuming and installation dependent. Complex arithmetic can 








be avoided by using a block diagonalizing transformation Uy 


rather than the diagonalizing transformation i feel, 
complex block diagonal matrix 


va fa[) 7} af 3] 


transforms A to a block diagonal real matrix 


1 ieee 


B= “AK =X “U “Ace 
ro} C28) Co} 28) 
= diag : yO * ss 
5 a is ea 


The matrix U&# is also real. If U is given as 


U= [(py+9,5): @y-494,5):-- +: tas): 74,5) 


then U, =U is 
Uy Sy ey oh ene a aed 


The eigenvalues are then given in terms of Uy as 


ee ae 
wou au wot), = Bay 
Ath eee 

[UL Ag UL Jig * 96 >i 


If the complex matrix ™ is represented as 


KH KH, i: 





The 


(A. 


(A. 


(A. 


(A. 


(A. 


(A. 


(A. 


18) 


19) 


20) 


21) 


22) 


23) 


24) 





if 


and #4 ~~ is represented as 


Pog = = 
HK = (J (J ; 

Ow), + (OD, (A.25) 
then the real and imaginary parts of the eigenvalue derivatives 
can be evaluated using real arithmetic as 


3 J -4 es Say 
go i" PM UL Ag OF dag 


a ced 
~ GO My Ag OF 14) a5 


+10 Uy Meby w+) + ow Ue ast, f+) ea: 


ries hat 


(A. 26) 


A.1.2 Differentiation of the Euler-Lagrange Equation Matrix 


An analytic formula for the derivative of the Euler- 
Lagrange equation matrix with respect to elements of its con- 
stituent matrix is quite complex due to the presence of several 
matrix triple products such as H'a H and Gp tc!. For this 
reason, the derivative of the Euler-Lagrange matrix is formed 
by a difference approximation to the derivative. The centered 


difference approximation to the first derivative is used. 


Piers £(xth) - £0ch) 


The resulting approximation is exact when used to evaluate the 
derivative of the Euler-Lagrange matrix with respect to any of 








its constituent matrices. The one exception is elements of 
the B_ matrix. 


A.2 LQR SENSITIVITY MINIMIZATION [56] 
A.2.1 Problem Statement 
Consider the system 
x = Fx + Gu+TIw 
(A. 27) 


y = Hx + Vv , 


where the matrices F, G, IY, and H contain parameters the 
values of which are uncertain. If these parameters are con- 
sidered as Gaussian random variables, a Gaussian random vector 
may be formed of which they are the components. This vector 


is specified 
E(z) = Zz, (A. 28) 


where the components of z, are the nominal values of the 


parameters, and by 
E 
E[(z - Zz) ( a, ] (A. 29) 


is a covariance i. “rix which is assumed known. Equation (A.27) 


may then be written as 
x = F(z)x + G(z)u + P(z)w (A. 30) 


HZ Sv y 


< 
tt 


x(0) = 0. 





The control is defined as 
u = -Cx (A. SDP) 


where the values of C may be left free or defined by function- 
al relationships to other parameters of the system. For this 
definition of u_ to be valid, Eq. (A.27) generally is required 
to describe and augmented system that includes the plant and 
the compensations. The matrices F, G, H, and [I then have 
to be defined accordingly. A quadratic PI for this system is 


2 
£ 
J = lim =a) B(x! ax + u' Bu) dt <2 (A. 32) 
13 
Ete £ 
) 


In this expression, the expected value is taken over the prob- 
ability distributions of x and u that are derived from both 
the distributions of the random process w and of the random 
vector z. Note that w is the process noise of the augmented 
system. Since w and z are independent, the expected value 
of a function of x and ui 4s 


E[f(x,u) ] f[x(w,z), u(w,z)] p[x(w,z), u(w,z) ]dxdu 


po 
Ww 
wv 
ad 


Ep (w) dw -f epcoe ‘ 


zZ 


: 
J 


where 


E. “f 8000 p(z)dz 


yA 





is the expected value over the distribution of z and 


EW = fen p(w) dw 


W 


is the expected value over the distribution of w. 


A free parameter vector @ is defined by the designer. 
This vector consists of q parameters of the system matrices 
that can be varied by the designer. 


The problem is now stated as follows. Given the system of 
Eqs. (A.30) and (A.31) in which the system matrices are functions 
of a variable parameter vector defined by Eqs. (A.28) and (A.29) 
determine the value of the free parameter vector 9 so that 


the PI of Eq. (A.32) is minimized. 


A.2.2 Derivation of Equations [56] 





The cost index Eq. (A.32) may be transformed to the form 


t 


! 
J=tr (A+C!BC) lim a [X_(t) + 8X(t) Jat 
oe te n 
O 


ts 
(A.34) 


= tr[(A+C Bc) (X, + 8X) ] 


where x and 6X are the time averages over all time of x 
and 6X, respectively, and A is the nominal state covariance 
matrix and 6X is the addition covariance due to parameter 
uncertainties. 

For a stable system, Xn and 6X tend to constant values 


i and 65X(~). Since the averaging in Eq. (A.34) is per- 


~ oe eee ee = 





formed over a large time interval, it can be assumed that 
xX * X (=) and 8X > &X(~). The PI therefore becomes 


J = 


tr [A+ CBC] [X, (@) + 8X(~) ] (A.35) 


where 
J, = tr[(A+C'BC)X, (*)] 


is the nominal PI, and 
Jy = tr[(A+C'BC)8X() ] 


is the additional PI due to the parameter variations. The PI 
of Eq. (A.35) can be minimized computationally by a two-step 
sequence: (a) the matrices Xt) and 6X(~) are found for 
a given value of the free parameter vector, 9; the PI that 
corresponds to this value is then determined; and (b) the 0 
vector is modified in a direction that decreases J. This 
sequence is repeated until the decrease in J in one cycle is 


less than a predetermined value. 


The two parts of this sequence are independent and computer 


programs for each one can be developed separately. 





APPENDIX B 


INTEGRAL CONTROL WITH FEEDFORWARDS - AN EXAMPLE 


B.1 STEADY-STATE SCHEDULING 


The philosophy of steady-state or reference point design 
can be illustrated with a simple linear example. Many problems 
with reference point selection in nonlinear, real systems are 
caused by the nonhomogeneity of the governing equations, rather 
than the nonlinearity, so a nonhomogeneous linear equation will 
show the important effects in an analytically tractable form. 


Consider Figure B.1: 


k 
f mk = f - kx 
yeax 
44 
Let m=ke=l1 
X¥ + x = 2u = gu 


Equilibria: Xog = 2Uc¢ (independent of any control law) 


Figure B.1 Example System 


—' 
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Steady-State Schedules 





There is a schedule to indicate how the output is to be 
programmed as a function of input quantities from the operator 


or the environment; this is shown in Figure B.2: 


a is an input specified by 
"pilot" or environment. By 
SS definition, there is some 

desired output from the system 


as a function of a. 


These three graphs show 
schematically the desired 


reference point of the system. 


Ss Each of the states and con- 
trols is really a function of 


only a _ but this graphical 


aN 


: representation is more conven- 
SS jent and computationally more 


efficient. 


¥Ss 


< 


Ss 


Figure B.2 Convenient Reference Schedules 











An alternative is to represent each of the states and 
controls as a function of a explicitly, for example, in 
Figure B.3: 


Seg Uss ; 
a a 


Figure B.3 Alternative Schedules 


But this would be much less computationally efficient. In 
complex systems, convenient, efficient and accurate reference 
point determination is a major design goal which can be achieved 
only through a physically meaningful approach and accurate sys- 
tem data. 


Returning to the example, a control system can be designed 
so that the output reaches a desired set point in a well 
behaved way. Consider the classical design using position and 
rate feedback (both assumed measured): 


Yes C(x - Xq) + Cv = 5) (B.1) 


Let 


Substituting into system equation: 


X + 2x + 2x - (xq) = 2(ug) (B.3) 








If the reference point is an equilibrium, 
cao 2u, 
% + 2x + 2x = 4uy 

and in steady-state, 


x = 2uq = Xq 


(B.4) 


(B.5) 


(B.6) 


The transient response to reference point changes and dis- 


turbances is well behaved in a classical sense (€ = 0.707). 


The behavior would be different if the actuator effectiveness, 


g, was specified incorrectly; the equations in terms of g 


explicitly are: 


x * gx * (2 +8) x - 5 Ba BOG the 


i} 
N 


Let g + € 
K+ (2ee)k + (2 + 5)x - (1 + 3)xy = (2+8)uy 
The schedule was based on the incorrect g _ so that 
ia = Zz Ug 
K+ (2+e)k + (2 + S)x = 2(2+e)uy 
and at equilibrium: 


. #{2%e) 


cc a 2 il | 


14 


(B.7) 


(B. 8) 


(B.9) 


(B.10) 


SR I nr ee are 





we ge aera . TT I 





and the steady-state "hang-off" is 


2€ 





Ax = Hog > Ma * aoe By (B.11) 
A r — = 
Au = Ugg Ug wie M4 (B.12) 


The transient response is also degraded but conservative 
control design should make this effect unimportant. It is 
possible that the steady-state offset could cause significant 
problems for high performance systems requiring accurate 
tracking. 


It can be shown that nonlinearities produce effects analo- 
gous to modeling errors so it is not possible to specify a 
reference point for a physical system without some inaccuracy. 
However, since physical laws can be applied to reference point 
schedules, a good design technique can reduce the error associ- 
ated with the approximation. 


B.2 INTEGRAL CONTROL FOR REFERENCE POINT TRIM AND LIMITS 


Suppose the control law in the previous example is modi- 
fied. Let 


Ue My - 5 (X-Xq4) -v- C,b (B.13) 
bax = Ky (B.14) 


Differentiating u and replacing in the system equation 





yields 
X° + 2% + 2x + 2C,(x-xg) = 0 (B.15) 
15 


— ——— 


At equilibriun, 


2. a XG (B.16) 


Seal 
1 di 
Au = (F 2 x) Xq (B.18) 


LE Cy is chosen so that the system is stable and the 
variation in g is not so large as to cause instability, then 
the state, x, will achieve its desired value independent of 
any system parameters. However, there still will be a hang-off 


in t. 


Many problems arise with this implementation. Most are 
due to the fact that transient behavior is degraded by the 
addition of a slow closed-loop root. Feedforward of the 
reference point is used in practical systems to achieve reason- 
able transient performance. Also, by raising the order by 
one, the parameter sensitivity has been increased. The stabil- 
ity region for variations in g and Cy is shown in Figure 
Buds 


The control law specified in Eq. (B.13) includes a control 
reference feedforward. However, it was shown that x 
achieves its desired value without any knowledge of the equi- 
librium u. Why should this feedforward be included? Consider 
the Laplace transform of the system with integral control: 


gC 


Ss 


gC, 
Ss 








(s? * gs * ] * . * 1) x(s) = +§ Xq4 + ug (B.19) 
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Figure B.4 Stability Limits of Integral System 


or for g¢ = 2, 


(s? +2574 2s 26) x(s) = (5 + 2G.) x, + Zsu (B.20) 
I 17d d 
With no feedforward, 
Ug = 0 
and (B. 21) 
x(s) _ Ss 2C, 
= 
Male) (gs? = a5” + 25 + 2C,) 
For a reasonably decoupled integral root (say 10x slower) 
x(s)_ 1 ( ie aaa (B. 22) 
xq(s) 5? ee a eal Sy 


Now, let the feedforward be modeled from the assumed 


equilibrium condition as follows: 


17 





Xd 
uy = = (B. 23) 
m 
and (assuming decoupling), 
2 + a. 2g. 
eS eae ae 
1 Ce eee, es aa a (B. 24) 
Mas) 4 ai ze © 2 tea 
Notice for g = 25 ‘that 
a 2 
po ena ts Col erate (B.25) 
Xq ts a $ 25-4 2 


and the transient response is independent of the integral con- 


trol for any input. This is a desirable property. It will 
not be possible to achieve this in general because of modeling 





inaccuracies and plant nonlinearities. However, it is intuitive- 
ly clear that the more accurately that Sa approximates g, 

the more the system behaves like the transient regulator design 
without the additional integration. 


The step response for various values of integral feedback 
gain is compared to that of the system with control reference 
feedforwards in Figure B.5. 


The moral of the story is that integral control provides 
a method for compensating small reference point specification 
errors and plant nonlinearity, but cannot replace carefully 
devised reference point specification algorithms because of 
degraded transient behavior. 
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Figure C.2 Engine Corrected Airflow Base Schedule as a Func- 
tion of Power Lever Angle and 
Engine Face Temperature 
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Figure C.3 Airflow Roll-off Factor Schedule for 
Low Speed, Low Power Modification 
of Scheduled Engine Airflow 
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Figure C.4 Airflow Corridor Required for Inlet 
Compatibility at All Power Conditions 
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Figure C.5 Fan Speed Schedule as a Function of 
Scheduled Engine Airflow 
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Figure C.6 Fan Match Point Schedule Consisting of a Base 
Schedule as a Function of Engine Airflow and an 
Additive Bias Shift for Supersonic Operation 
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Figure C.7 Engine Pressure Ratio Schedule as a Function 
of Scheduled Fan Speed 
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Figure C.8 Scheduled High Rotor Speed as a Function of Power 
Lever and Engine Face Total Temperature 
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Figure C.12 Nozzle Area Scheduled as a Function of 
Power Lever Angle and Engine Face Total 
Temperature 
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Figure C.13 Rear Compressor Variable Vane Angle Schedule as 
a Function of Corrected Compressor Speed and Flutter 
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Figure C.14 Compressor Inlet Variable Vane Schedule as a 
Function of Corrected Fan Speed with 
Flutter Boundary on Maximum Deflection 
as a Function of Engine Face Total Pressure 
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APPENDIX D 
LINEAR MODELS 


D.1 DETAILED LINEAR MODEL ANALYSIS 


Linear models were generated at a group of operating condi- 
tions throughout the flight envelope [1]. Operating conditions 
were specified sequentially in groups. Operating point speci- 
fication was determined by insight into the behavior resulting 
from the dynamic analysis of the previous groups of models. 

The 42 operating points chosen are shown in Table D.1. The 
analysis of the dynamic behavior coupled with the detailed 
steady-state analysis contained in Chapter V produced the choice 
of design and evaluation points for regulator synthesis. The 
details of the linear model analysis of the matrices, generated 
by Pratt §& Whitney Aircraft and included in their report [1], 
are presented below. 


Sea Level Static Intermediate Power Models 


Linear models were generated at a series of operating 
points at sea level static conditions. For intermediate power, 
the root locations of the model are shown in Table D.2. It is 
possible to associate a state variable with each root by inspec- 
tion of the eigenvectors of the system. A correspondence can 
then be found to the models of pressure, temperature and flow 
time constants included in the nonlinear simulation. In this 
way, the structure of the response can be determined in terms 


of predominant response variables. 


The modeled state response to initial conditions for five 
of the engine variables has been calculated to illustrate 
the correspondence of eigenvector components and modal responses 
in the linear system. In Figure D.1, the time history of the 


Table D.1 
Operating Points Chosen for Dynamic/Static Analysis 
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five states is calculated for unit initial conditions in each 
state. In the first column, fan speed is initially one unit 
away from the trim speed. As the system returns to equilibrium, 
it is observed that the augmentor pressure lags the fan speed. 
The compressor speed is nearly unaffected. The increase in 

fan discharge pressure is transmitted into the burner because 
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Table D.2 
Eigenvector Association of Dynamic Roots of 16th Order 
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Figure D.1 Indicial Response of Fifth Order Linear Model, 
Intermediate Power, Sea Level Static 
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of the nearly constant compressor pressure ratio. As the fan 

speed moves to equilibrium, the burner pressure follows. For 
an initial condition in compressor speed, the response is 
nearly a simple exponential decay. The fan speed and augmentor 
pressure move in phase in response to the compressor speed. 
Burner pressure follows the compressor speed decay. For a 
temperature and burner pressure initial condition, the decay 

to equilibrium does not excite any motion in the fan, compres- 


sor speed and augmentor pressure states. 


The dynamic character of the linear system is evident 
from the eigensystem analysis and the indicial response plots. 
The eigensystem shows a root with a 7 msec time constant. This 
root represents the flow equilibration in the burner as illus- 
trated in the 5th column of Figure D.1. Similarly, the gas 
temperature response of the system at the turbine inlet also 
has a fast time constant relative to the rotor speed roots. 


The complex root near |s| = 5 rad/sec has response charac- 





teristics involving fan speed and an out-of-phase component. 
The real root is associated with the compressor speed 

state. The behavior predicted by the eigenvectors is 

directly observed in the linear simulation as would be 
expected. A block diagram of the linear system can be con- 
structed showing the important interactions. Figure D.2 shows 


this diagram. 


The correspondence of the linear models to the engine 
response generated by the nonlinear simulation has not been 
thoroughly testede The dynamic response of pressures, 
temperatures, and rotor speeds predicted from the linear analy- 
sis, in general, seems adequate. However, the strong dynamic 
interaction between the augmentor pressure loading on the fan 
disc and the flow equilibration in the duct/augmentor volume 
is certainly not predicted from simplified physical analysis. 
The adverse effects of modeling inaccuracies are greatly re- 
duced within a closed-loop feedback system such as the linear 
quadratic regulator. Thus, precise models in closed-loop 





Figure D.2 Linear System Block Diagram 


design are not necessary to achieve excellent performance as 
long as the synthesis procedure is undertaken with a recogni- 
tion of the approximate nature of the representation. The 
linear models are used in LQR design in this spirit and with 
the firm understanding that they are merely a means to an end, 
namely, a practical controller for a real engine, 


The analysis can be used on the output response of the 
system. Figure D.3 shows the indicial output response (i.e., 
the response to unit initial conditions in each state). The 
primary components of the thrust are the fan speed, compressor 


speed, and augmentor pressure. 


Extension of Results to the Sea Level Static Operating 


Line 





The applicability of the linear model to variation in 
power level is the next issue. In Figure D.4, the frequency 
response at idle and intermediate power is shown for the fuel 
flow to low rotor speed transfer function. At intermediate 
power, the low rotor is driven by the fuel flow through the 
cascaded high rotor speed lag. As a result, the transfer func- 
tion is nearly first order with a response showing only the 


high rotor dynamics. 
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Figure D.3 Output Indicial Response 


At idle, the response of low rotor speed to fuel flow 
shows the effects of the compressor lag. This indicates that 
a mode of independent low and high rotor motion, suppressed 
at 1igh power, can be excited at low power. Also, the high 
rotor root has only one-third the response time at idle. 


An intuitive explanation of this behavior is suggested 
from engine operation. At high power, the fan turbine operates 
at a constant flow parameter. Changes in high turbine inlet 
temperature are reflected as excess work at both turbines 
instantly (or nearly so). Fuel energy is converted to torque 
for both spools and the response of the system to fuel flow 
is a characteristic lag of the combined rotor system acceler- 
ating (i.e., the two spools act as a single rotor system). 
Pressure disturbances in the augmentor appear only at the com- 
pressor inlet because the choked turbine eliminates any back 
pressure sensitivity. 
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At low power, the fan turbine is operating in an unchoked 
condition. Increases in turbine inlet temperature cause 1in- 
creases in the inlet flow parameter and torque. The response 
of the core engine is determined by the high spool alone. The 


reduced effectiveness of converting temperature to torque 1n 


the fan turbine causes the response to be slower. At idle, 
the exhaust nozzle is unchoked. Thus, pressure disturbances 

in the augmentor dissipate quickly through the nozzle discharge 
with a significant increase in flow equilibration rate. At 
low power, the engine is quite nonlinear because of the many 
types of aero-thermodynamic couplings which occur at low flow 
conditions. A linear model determined from an offset derivative 
algorithm will be step size dependent. The control design at 
these conditions must operate satisfactorily in the presence 

of the extremely variable behavior and handle transitions to 
regimes where the behavior is more amenable to a linear differ- 


ential model. 


Extension to Subsonic Operating Line Models 





The linear response of the engine at sea level static power 
conditions is characterized by a simple linear response at 
conditions from above idle to intermediate power. Near idle, 
unchoking of the nozzle and low turbine cause nonlinear behavior 
characterized by a slower, uncoupled response of the two rotor 
speeds. The modes of operation at various power conditions and 
altitude points are related to this behavior. The models can 
be parameterized with an ambient variable when a correspondence 


is found. 


Figure D.5 shows the engine operating envelope. Transi- 
ent operation is limited to the left by low pressure constraints 
and minimum airflow requirements. To the right, operation is 
restricted by maximum pressure, temperature and airflow require- 
results in a narrow corridor in which transitions 


ments. This 
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in power condition are permitted. In these regions, the dy- 
namics of the plant must be understood so that practical and 


minimally complex controls can be implemented. 


From sea level static, lower points in the transient cor- 
ridor correspond to altitude, subsonic conditions. Higher in 
the corridor corresponds to higher Mach number and low altitude 
points. Two series of operating line models were generated 
at subsonic altitude points. At 30K/0.9, the ambient pressure 
is roughly one-half of sea level static and at 45K/0.9, it is 
roughly one quarter. Higher than 45,000 ft, the subsonic 
operating envelope is limited by pressure constraints to in- 


clude only intermediate power. 


Part power models for 45K/0.9 and 30K/0.9 are shown in 
Table D.3. At both conditions, the nozzle and fan turbine 
are choked for the series of power conditions. At high power, 
the pressure and rotor speed roots are uncoupled. Since the 
rotor speed time constant at a fixed, nondimensional match 
point is proportional to 68V@ [69] and the pressure equilibra- 
tion time constant is approximately proportioned to /Y6 [70], 
the coupling between these states should be weaker at higher 
altitudes and low Mach numbers. For the high rotor speed 
time constant, these relations are accurate. For the low rotor 
and pressure roots, the coupling terms in the linear models 
appear to be important, particularly at low airflows. The 
root locations reflect this coupling. Thus, at high power, 
the system has two uncoupled dynamic roots corresponding to 
the rotor speed response. This behavior is reflected in the 
30K/0.9, 45K/0.9 and 65K/0.9 models. 


At lower powers, a different behavior becomes apparent 
at 30K/0.9. In this case, there is a complex root associated 
with the spool response. The coupling is probably due to the 
load on the compressor face caused by changes in fan speed 
and discharge pressure. This effect is not dynamically im- 
portant at higher pressures Or airflows because of the decreas- 


ed sensitivity of compressor speed to compressor entrance 
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perturbations. It should be noted that the dynamics of the 
States are closely coupled at these points. This contrasts 
with the response of the system at the higher power conditions 


where the system is more uncoupled. 


The conclusions concerning altitude points can be sum- 
marized as follows. The high power behavior predominantly 
corresponds to uncoupled rotor speed dynamics. At low power 
conditions, moderate altitudes and Mach numbers, the spool 
responses couple to form a moderately damped complex pair. 

This behavior is probably produced by the increased sensitivity 


of the high pressure spool to compressor face disturbances. 


Several other subsonic points were investigated. These 
points were generally low altitude, low Mach number at inter- 
mediate and idle power. Figure D).5 shows the low turbine 
Choking parameter Tr4.5/Tr2 as a function of airflow. The low 
turbine is observed to run unchoked in a region of the flight 
envelope at low altitude and Mach number conditions. The 
eigenvalues of the models generated to investigate this regime 
are listed in Table D.4. These models demonstrate the response 
predicted by the unchoked turbine behavior. It can be seen that 
the character of the roots changes significantly, even at points 
which are quite close together in the envelope (e.g., 30K/0.9, 
30K/0.6 and 20K/0.3). This is indicative of the sensitivity 
of the linear dynamics to the chosen operating point; or, 
stated another way, the nonlinearity of the plant in this 


region. 


Extension to Transonic, Low Altitude Points 





The behavior of the engine at conditions corresponding to 
higher inlet pressures and temperatures in the upper half of 
the transient corridor was investigated with models along the 
operating line at 10K/0.9 and intermediate power at OK/1.2. 


These results are shown in Table D.5. 
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Table D.4 


Eigenvalues of Miscellaneous Subsonic Linear Models 


83/20K/0.3 24/20K/0.3 20/10K/0.6 24/30K/0.6 





* 
Dominant roots 


The behavior of the models at 10K/0.9 is characteristic 
of the high power system at sea level static conditions. The 
results also indicate that the model does not change along the 
operating line. This condition operates at higher corrected 
airflow at idle than lower altitude, lower Mach number points 
and both the nozzle and turbines are choked. 


The high ram pressure condition, 0K/1.2, has a limited 
power level range because of airflow constraints near inter- 
mediate power. The response of the engine is similar to the 
sea level static, intermediate behavior. The response times 
of the roots are about 30% less than at sea level static. 
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Table D.5S 


Eigenvalues of Transonic Operating Points 


83/10K/0.9 | 67/10K/0.9 | 52/10K/0.9 | 36/10K/0.9 | 20/10K/0.9 | 83/0/1.2 


-21 + .8j | 
-12 | 
-5.6+1.7j* 
<1 1% 
=? 
-0.8 





* 
Dominant roots 


The conclusions concerning the behavior in this portion 
of the envelope are summarized below. The eigensystem of the 
models in this region of the envelope is similar to the sea 
level static, intermediate power. The response times of the 
engine decrease as a function of the position in the envelope 
with the sea level ram point being about 30% faster than the 
sea level static response. 


Extension to Supersonic Conditions 


Several supersonic points were chosen to investigate engine 


dynamic response. These points range from very low to very 
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high ram pressures and at higher than standard inlet tempera- 
tures. Airflow limits constrain the engine to intermediate 
power (and above). Thus, these points lie to the right of the 
transient Operating corridor. The control requirements at 
these points are compensation of disturbances due to augmentor 
lights and inlet distortion. 


Three groups of models were chosen. These were low pres- 
sure, mid-pressure, and high pressure points in the right side 
of the envelope. The results from the model analysis for each 
group is presented below. 


The low pressure points consisted of 75K/1.8 and 58.5K/ 
Z-15. The model at 75K/1.8 is very similar to that at 65K/0.9 
and other subsonic altitude points. The rotor speed responses 
are coupled. The pressure lag of the augmentor is considerably 
faster. Control of pressure pulses at this condition due to 
augmentor lighting can be detected using a pressure measurement. 
A disturbance in the augmentor propagates quickly to the fan 
exit station, well before the spool speeds have responded. 

At 58.5K/2.15, the operating point has the same density as 
30K/0.9 but at a higher pressure and temperature. The correct- 
ed airflow at this point corresponds to approximately a PLA 

of 35° at 30K/0.9. The 36°/30K/0.9 model is similar to the 
model at 83°/58.5K/2.15. 


The mid-power supersonic points were 65K/2.5 and 40K/1.8. 
At 40K/1.8, the model has coupled rotor roots. This point has 
roughly the same pressure as sea level static and roughly the 
same temperature as the O0K/1.2 point. The behavior is quite 
different than either of these points. It lies along a con- 
stant density line corresponding to 20K/0.9 conditions at a 
PLA of approximately 40°. No models were generated at this 
point but extrapolation of the low power, 30K/0.9 behavior 
would indicate a system response of the type exhibited by the 
40K/1.8 model. At 2.5/65K, the operating point has the same 
density as 30K/0.9 at extremely high inlet temperatures. The 
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airflow is low, corresponding to 28°/30K/0.9. The model does 
not exhibit this type of behavior; rather, it is more typical 


of the sea level static response at low power. 


The high pressure points consisted of 20K/1.8 and 40K/2.2 
at intermediate power. The 20K/1.8 corresponds to the density 
at OK/1.2 with higher pressures and temperatures. The behavior 
of these models is nearly identical. The 40K/2.2 operating 
point has the same pressure as OK/1.2, and roughly the same 
temperature as 58.5K/2.15 at a lower density than 0K/1.2. 

The response has a coupled rotor system. This behavior does 
not correspond to any other high pressure system. It has the 
same Characteristics, however, as 58.5K/2.15 linear system. 


The conclusions important in the analysis of the super- 
sonic operating conditions are as follows. The low pressure 
supersonic points are similar to the corresponding constant 
density operating points at subsonic altitude points with the 
same corrected airflow. Higher temperatures tend to destroy 
this correspondence. 


Supplemental Models 
A linear model was generated at the 45K/0.9 flight point 


at maximum (full augmented) power. This point is representa- 
tive of the flight conditions where pressure spike suppression 
during augmentor ignition is important. The low ambient pres- 
sure and fan stability margin at intermediate power make this 
flight point critical to satisfactory augmentor operation. 

A comparison of the intermediate and maximum power model roots 
are shown in Table D.6. During augmentor operation, the engine 
airflow is nearly constant and the operation of the core engine 
and fan should be the same throughout augmented operation. 

It is important to analyze the differences in the transient 


behavior to assess the impact of differences on performance. 
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Table D.6 


Eigenvalues of Intermediate Power and Maximum 
Power Systems at 45,000 Ft, 
Mach Number = 0.9 


-19 
-6.3* 
-1.2+0.4j* 





* 
Dominant roots 


Recalling the intermediate model at these flight condi- 
tions contained an uncoupled set of real roots, the fully aug- 
mented roots exhibited lower frequency, coupled, spool response. 
Since augmentor control is primarily pressure modulation, the 
pressure dynamics of both mdels are similar. 


A linear model was generated at sea level static idle 
conditions with the full specified level of componert deterior- 
ation, horsepower extraction and bleed flow. It was determined 
that sea level static idle conditions were the most sensitive 
to installation and deterioration effects. The results are 
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shown in Table D.7. Clearly, the impact on the engine dynamics 
is small while the change of the control effectiveness is more 


pronounced. 


Table D.7 


Eigenvalues of Models at Sea Level Static, Idle Power 
With and Without Full Installation and 
Deterioration Effects 


20/0K/0 20/0K/0 
NONE NOMINAL 
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D.2 SELECTED LINEAR MODELS 


5x = FSx + GSu 


H8x + Du 


by = 
Note 1: For state, control and output ordering, see 
printout and symbol table. 
Note 2: Models are normalized by the baseline operating 


point (with the exceptions shown below) so that 


units of all quantities are dimensionless. 
(DP25QT) ony = 1.2/0 


(DP25SC) voy = 71.10 


Note 3: Table D.8 lists the linear models included. 











a 


a * 
- 
Re 


Table D.8 





Index of Linear Model Data 
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D.3 SELECTED REDUCED ORDER MODELS 
D.3.1 Definition of Reduction Matrices (see Chapter IV) 
Modal Decomposition 
5x = F8x + Gdu 
5x = Tz 
FT = TA 
where 


A is block diagonal 


Partition: 


= ' ie = 
24 Ay 0 24 =) 
24 0 A, Z4 E, 
5x) Tay 1 a9 1% 
5x4 Ta, ! T22 JL22 
5x, 
By = [H, H,]|----] + Déu 
5x, 
148 











Reduction: 


6X, = F 5x 


1 + G,5u 


1 


dy = H).6x, Bs D Su 


= Ys * 
5x, H 5x, + D*5u 


Table D.9 shows the equivalence of the above matrix names 
and the listing. 


D.3.2 Normalization 

Matrices are normalized as in Section D.2. The eigensystem 
shown in not normalized and matrices listed at the end are the 
final reduced, unnormalized system. 


D.3.3 Case Index 


Table D.10 lists the cases included. 
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Table D.9 


Index of Reduction Matrices 
(In order of appearance in printout) 


LAMBDA1 


Th 


T 
LAMBDA2 


T,, INVERSE 

FT 

LAMBDA2 INVERSE 
GR 

H* 

HR 

D* 

OR 


Unnormalized Output 





a eee SF 











Table D.10 
Index to Model Reductions 


PLA ALTITUDE 
CASE (DEG) (FT) MACH NO. 
83 0.9 


45,000 
83 30,000 0.3 
0 
0 
30,000 
0 
0 
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EXAMPLE OF TRANSFER FUNCTION DATA ON 16TH ORDER LINEAR 
MODEL 


PLA = 52° 
ALTITUDE = 
MACH NUMBER = 0 


ELEMENT UNITS VARIABLE SYMBOL 


Fan Speed 
Compressor Speed 


Compressor Discharge Pressure 


Interturbine Pressure 
Augmentor Pressure 

Fan Inside Temperature 

Duct Temperature 

Compressor Discharge Temperature 
Burner Exit Fast Temperature 
Burner Exit Slow Temperature 
Burner Exit Total Temperature 
Fan Inlet Fast Temperature 
Fan Inlet Slow Temperature 
Fan Turbine Exit Temperature 
Duct Exit Temperature 
Augmentor Temperature 


Fuel Flow 
Nozzle Area 
CIVV 


0° y aT ee on. 2 BE RT 
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APPENDIX E 


LINEAR QUADRATIC REGULATOR DESIGN DATA 
(OPTSYS PROGRAM) 





Plant: x = Fx + Gu 

Outputs: y = HY x = DY u 

Performance l oe + 

Index: J = a | [x Ax + u Bu] dt 
(e) 


Low Rotor Speed 

High Rotor Speed RPM 
Augmentor Pressure psia 
Fuel Flow Tbm/hr 
Burner Pressure psia 


Fuel Flow Command lbm/hr 

Nozzle Area Ft? 

CIVV deg 

RCVV deg 

Bleed Flow Fraction Compressor Flow 


Thrust lbf 
Airflow 1bm/sec 
TT4 “AN 

Fan Surge Margin Per cent 
Compressor Surge Margin cent 
Ap/p Parameter 


Ap/p Average 





an eo 
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LQR BASELINE DESIGNS 
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E.2 SUPERSONIC EVALUATION POINT DATA 


This section contains closed-loop unnormalized eigensystem 
using off-design regulator gains from gain schedule algorithm 


(for open-loop data see Appendix D). 


ALT fe 
CASE — DESCRIPTION 


I 58.5K Middle Pressure 
II Low Pressure 
III High Pressure 


(16th order linear model plus one state for fuel flow lag) 
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APPENDIX F 


SENSITIVITY DATA FOR LINEAR QUADRATIC REGULATOR GAINS 
AND GAIN SCHEDULES FOR MULTIVARIABLE CONTROLLER 


F.1 UTILIZATION OF SENSITIVITY DATA TO SIMPLIFY CONTROL LAW 


The controller gains resulting from LQR designs at the 
six design points (Appendix E) are compiled in Table F.1. 
These gains were smoothed and simplified using the sensitivity 
data in Section F.2 to form the final set of implemented 
control gains shown in Table F.2. 


F.2 SENSITIVITY DATA 


Sensitivities of the closed-loop eigenvalues are given 
as a function of each element of the regulator gain matrix and 
each diagonal element of the state and control weightings. 


to parameter, 6. 


The sensitivity of the eigenvalue, A). j? 


les 
is defined as follows: 


where 8, is the baseline value of the C, A, or B 


matrices. 
The sensitivity S5 is a complex vector which is associ- 


ated with a particular wigerralue r;- The printout is organ- 
ized to show all the sensitivities of a particular matrix 
element of the C, A, or B matrices in a group. The 
Baseline Parameter Value is the value, 85° in the matrix. 
The Pole Location column contains values of the five eigen- 
values, Ass in the fifth order linear design mode. The 
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Table F.3 
Gain Component Symbol Index to Tables F.2 and F.3 


Main Burner Fuel Flow Tbm/hr 
Low Rotor Speed rpm 
High Rotor Speed rpm 
Augmentor Entrance Total Pressure psia 
Requested Fuel Flow 1bm/hr 
Burner Pressure psia 
Jet Area ft? 


Compressor Inlet Variable Vane Angle deg 


Rear Compressor Variable Stator Angle deg 
Bleed Percentage Flow fraction 





Table F.4 


Index to Sensitivity Printout 
(Design Point Data) 


PLA ALTITUDE 
CASE (DEG) (FT) MACH NO. PAGE 
83 0.9 314 


45,000 
83 30,000 
83 0 
83 0 
20 30,000 
20 0 
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Sensitivity data is organized into two columns containing the 


Real and Imaginary part of the complex vector: 


i Oe 4 
S ° @ SB eee 
6. fo) 00; 
Also, the Magnitude and Phase of Ca is given. Note, if 
J 


8, = 0, then 8, = 1 is assumed. 
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APPENDIX G 


INTEGRAL CONTROL DESIGN DATA AND EXAMPLE DESIGN 


G.1 DATA FOR INTEGRAL CONTROL DESIGN 


The integral control designs were carried out at the six 
control design points according to the procedure described 
in Sections 4.6 and 5.8.3 of Volume I of this report. The 
data used in the designs is included in this section. The 
matrices used in the general development presented in Section 
4.6 are shown here for the conditions listed in Table G.1. 
Table G.2 reviews the general equations to define the matrix 
outputs presented. A specific example is shown in Section 
G.2. The order of states, outputs and controls is given in 
Table 6.3. 


Table G.1 
Index of Steady-State Gain Matrices 


DESIGN PLA ALTITUDE | MACH NO. PAGE 
POINT (DEG) (FT) (-) 
83 45K 0.9 392 


83 30K : 394 





























Table G.2 
Integral Design Relationships 


F5x + GSu 


HBx + DSy 


0, 


Cy 8x + Su 





-(F+6C)7! Geu- 

H6u~ 
-1 

-(H+ DC.) (F + GC.) G+D 


Integral System 


BYs = Y ~ Ygp 





Table G.3 
Order of States, Outputs, and Controls 


SYMBOL UNITS 


Ny RPM 
No RPM 
PTOM PSI 
WF LBM/HR 
PT4 PSI 


CONTROL 


LBM/HR 
FT? 
DEG 
DEG 
FRACT. COMPRESSOR 
FLOW 


OUTPUT 


FN LBF 
WCF LBM/SEC 
TT4 








SMAF 
SMHC 
(AP/P) eyp 
(4P/P) ave 





*Not included in some models 
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DESIGN CYCLE WITH VERIFICATION SHOWING BIGEN- 


G.2 EXAMPLE OF 
EAR SYSTEMS 


VALUES OF AUGMENTED LIN 


Section 5.960 (volume t} 35 


discussed in 
point, intermediat© 


troller for one design 
(83/0/09). The design pro- 


e an acceptable transient 
of closed- 


The procedure 


used to design 4 con 


power, sea level static conditions 


ceeds by using LQR methods t° produc 
yr design (step Wes Next, the sonsitivi tte? 


reguilato 
ea used to eliminate gain 


Loop response ar 
not important in the controlled response 


integral trim designs 4 
are augmented to the original system. 
for the simplifie 


possible trim modes 


d transient regu a each of three “ 


are calculated for 
4, and 5- This £1n@ 
Linear description fo) 


system in Steps ie 1 check assures accept~ 
able stability with 
augmented plant. Ta 
tabular data representing 
culated by the software: Figures 
f the closed- Loop Yr 


an accurate 
ple G.4 lists the steps an 
the actual eigenvalue loc 
G.1 and G.2 graphically show 


oots as the design and evalu- 


Table G.4 












g3/0K/0 
by eight 


Baseline LOR synthesis at 


Regulator matrix reduced 
terms 

Augmented system with five inte- 
gral roots - 1 Dp/P trim mode 
Fill = Ap/P trim mode 

m7 Ap/P trim mode 
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Figure G.1 





© OPEN LOOP 
[[] FULL LOR - STEP 1 


© REDUCED LOR - STEP 2 





Open-Loop Poles, Full LQR Poles and Simplified 
LQR Poles (Step 1 and Step 2) 
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APPENDIX H 
TRANSITION CONTROL DATA 


The development of the transition control is discussed 
in detail in Section 4.5 (Volume I). The fifth order design 
models were used to calculate the numerical quantities used 
for the transition control design. This data is presented in 
this section. A critical step in the design is the partition 
of the dynamical state equations to determine the primary set, 
5x) and the secondary set 5x5. Since the RCVV's, CIVV's, 
and bleeds have only a small range of motion, these controls 
were not used during large accelerations in the explicit 
transition algorithm. Fuel flow and jet area were the two 
controls used for the acceleration control The fuel 
flow is the primary acceleration regulating variable. Jet 
area has a strong effect on engine response, even though its 
excursions are limited, essentially, to increases about the 
steady-state values. Thus, the physical character of the F100 
behavior was the determining factor in the choice of modulated 


controls. 


The primary state rates chosen were the fan speed and 
compressor speed. These were used because the continuous 
state asymptotic expansion will provide continuous reference 
trajectories for these values which is close to being phsyical- 
ly realizable. Also, the effects of the two controls on the 
speeds were independent as required by an invertible G, 
matrix. This fact is directly calculated from the equation, 
or physically justified by realizing that fuel affects the 
rotor speeds collectively while nozzle area affects the speed 
differentially. 


The definitions of the partitions of the fifth order 
design models are given in Table H.1. The specific values 
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ae 
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Table H.1 
Variable Index 


5x1 Fay | Fa 
5x, Pou» Bag 
wee eee” 
5x 

by = [H, | 
<a 5x 
HY 2 


Low Rotor Speed 
High Rotor Speed 


5X5 


Augmentor Pressure 
Fuel Flow 
Burner Pressure 


* 
Not included in some 


psia 
Tbm/hr 
psia 





models 


ee a ooo 


5x, Gy 
----]| + |---|] 8u 
5x, G, 
G 
+ Dj du 
DY 


Fuel Flow Command 
Nozzle Area 
dy 
Thrust lbf 


Airflow lbm/sec 


Turbine Temperature 


Fan Stability Margin 
Compressor Stability Margin 


* 
Ap/p 


* 
(AP/P) avg 


— 


used in the design of the rate schedules are discussed in 
Section 5.9.2 and listed in Table 5.9 (Volume I). The for- 
mulas for the design of the rates are derived in Section 4.5 
and the formulas for data provided in the printout are shown 
in Table H.2. 


The data is presented for the equations specified in 
Table H.2 for the conditions used to derive the transient 
rates for the controller presented in Table 5.9. Eleven cases 
are presented as reference data and these cases are indexed in 
Table H.3. 


A diagram of the rate limited structure of the transi- 
tion model is shown in Figure H.1. For various size inputs, 
the normalized output is shown in Figure H.2. It is obvious 
that for small inputs, the response is determined by the 
regulator/trim control. For large transients, the transition 
control has authority until the steady-state condition is 
approached. The structure used is always globally state and 
does not use switching or other undesirable high gain nonlin- 


earities which can degrade overall stability. 
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Table H.2 
Matrix Index 


Intermediate Matrices: 


a 
Bay. * Soi Fig 


= 


F 1 


22 ~ &26) Fiz 


State Rate Compatability Matrix: 


State Initial Jump Matrix:* 


__ pl 
Sx,(0) = F, 


= 
[M+ G,G, ] as 


Control Initial Jump Matrix:* 


awd -1 -1 
Su(0) =G)° [I+F, F5°(M+G,G)")] a, 


Control Rate Matrix:* 


Spy os ' 
H, = Gy TFiy > Fygtl Oy 


Output Jump Matrix: 


=i 


4 =-1 -1 
= -H4F, [M+G,G, ] + DG, 


Output Rate Matrix: 


-4 <f 
1 ~ DG Fy, + (Hz + DG) Fy 2)M 





*Matrix included in data is underlined 


426 








Table H.3 
Index of Transition Data 


POWER LEVER 
ANGLE 
(DEG) 


ALTITUDE MACH 
(FT) NO. 


0 
30,000 
30,000 
45,000 
45 ,000 


oO On DD OT FP Won re 





Figure H.1 Basic Structure of Transition Control Logic 
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APPENDIX I 


EFFECT OF COMPONENT DETERIORATION, POWER EXTRACTION 
AND BLEED FLOW ON F100 TURBOFAN CONTROL 
SYSTEM DESIGNS 


Table I.1 
Index of Data 


FIGURE COMMENTS 


Altitude Points - Intermediate 
Power Data for Nominal Deter- 
joration and Bleed Flow 


Sea Level Static Performance* 
30K/0.9 Performance* 
10K/0.9 Performance* 





* 
Data for each effect is given as a function 
of power level 
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ALTITUDE 
(FT) 


OK 10K 20K 30K 45K 58. 5K 65K 


PRODUCTION CONTROL 
L LE 





Figure I.1 Per Cent Variation in Intermediate Power Net Thrust 
Due to Component Deterioration at Various Altitudes 
and Velocities from Digital Evaluation 


ALTITUDE 
(FT) 


OK 10K 20K 30K 45K 58. 5K 65K 


PRODUCTION CONTROL 
L L 





Figure I.2 Per Cent Variation in Intermediate Power Net Thrust 
Specific Fuel Consumption Due to Component Deterioration 
at Various Altitudes and Velocities from Digital Evaluation 


482 





—— = — = a 
ALGO LO AE CO tO . a a 





ALTITUDE 
(FT) 


OK 10K 20K 30K 45K 58. 5K 65K 


PRODUCTION CONTROL 
MULTIVARTABLE CONTROL 





Figure I.3 Per Cent Variation in Intermediate Power Net Thrust 
Due to Customer Bleed Flow at Various Altitudes 
and Velocities from Digital Simulation 
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THRUST 
(% CHANGE) 


THRUST SPECIFIC FUEL CONSUMPTION 
(% CHANGE) 






‘ 
on 


. VAR. CONT. - 5 NOM. DET. 
-10 . VAR. CONT. - NOM. DET. 


. CONT. - 5s NOM. DET. 
. CONT. - NOM. DET. 


-15 
20 30 40 50 60 70 80 90 


POWER LEVER ANGLE 
(DEC) 


. WAR. CONT. - 4 NOM. OET. 
. VAR. CONT. - NOM. DET. 

. CONT. = 4 NOM. DET. 

. CONT. = NOM. DET. 





20 30 40 50 60 70 80 90 
POWER LEVER ANGLE 
(DEG) 


Figure 1.4 Per Cent Variation in Steady-State Performance at 





Sea Level Static Conditions 
from Digital Simulation 
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THRUST 
(% CHANGE) 


THRUST SPECIFIC FUEL CONSUMPTION 


(% CHANGE) 











. VAR. CONT. - '; NOM. HPX 
. VAR. CONT. - NOM. HPX 

. CONT. - 4 NOM. HPX 
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Figure I.6 Per Cent Variation in Steady-State Performance 





. WAR. CONT. - 4 NOM. HPX 
. VAR. CONT. - NOM. HPX 

. CONT. - 5 NOM. HPX 

. CONT. - NOM. HPX 


THRUST 
(% CHANGE) 





20 30 40 50 60 70 80 90 


POWER LEVER ANGLE 
(DEG) 


0.4 


0.3 






(% CHANGE ) 


o 
nm 


THRUST SPECIFIC FUEL CONSUMPTION 


O.1 . VAR. CONT. - +; NOM. HPX 


. VAR. CONT. - NOM. HPX 
. CONT. - 5; NOM. HPX 
. CONT. - NOM. HPX 


20 30 40 50 60 70 80 90 


POWER LEVER ANGLE 
(DEG) 


Figure I.6 (Continued) 





10 






THRUST 
(% CHANGE) 


1 
Lawd 
o 


. VAR. CONT. - 4 NOM. BLEED 
-30 . VAR. CONT. - NOM. BLEED 


. CONT. - 4 NOM. BLEED 
. CONT. - NOM. BLEED 


-40 
50 60 


POWER LEVER ANGLE 
(DEG) 





THRUST SPECIFIC FUEL CONSUMPTION 


. VAR. CONT. - 4 NOM. BLEED 
. VAR. CONT. - NOM. BLEED 

- CONT. - 4 NOM. BLEED 

. CONT. - NOM. BLEED 





0 
20 30 40 50 60 70 80 90 
POWER LEVER ANGLE 
(DEG) 


Figure 1.6 (Continued) 


493 


nT ene . — nn RR ER I 











THRUST 
(% CHANGE) 


20 30 40 50 


© MULTIVARIABLE CONTROL 
——— PRODUCTION CONTROL 





60 70 80 


POWER LEVER ANGLE 
(DEG) 


THRUST SPECIFIC FUEL CONSUMPTION 
(% CHANGE) 


© MULTIVARIABLE CONTROL 
——— PRODUCTION CONTROL 


5 


60 





POWER LEVER ANGLE 
(0€G) 


Figure 1.6 


(Continued) 


494 


90 





APPENDIX J 


EXAMPLE OF LINEAR SPLINE SCHEDULE GENERATION 


Section 4.4 of Volume I describes the utilization of splines 
to fit operating data. This procedure is useful in generating 
gain and reference point schedules. An example of this proced- 
ure is given below for the generation of an optimum fan operat- 
ing line schedule. The program flow diagram is shown in Figure 
J.1. Schedule data for two dimensional tables are read in and 
weighted according to the importance of the operating point. 
Table J.1 shows a list of 62 steady state operating points 
represented by values of corrected fan speed and airflow. 
Higher weights are attached to low altitude points. An initial 
set of breakpoints are chosen and the data is partitioned into 
regions as shown in Figure J.2. The least squares problem is 
shown below: 


yy x1 1 0 0 ! 0 ay 
Yp *D> 1 0 0 0 by 
Yor 0 0 Xptl a 0 ao 
eho tJ.) 
Yq 0 0 scr 1 0 by 
0 , 0 0 : 0 Saad 


The parameter constraints representing continuity of splines 
at the breakpoints is represented as follows: 


Le = 0 (J.2) 
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Figure J.1 Spline Table Generation Program 
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Figure J.2 Spline Schedule Generated from Data 
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where 
= [a,, by» sees acs be] (I-35) 


and 


xX) 1 -X, =i 0 0 0 0 0 0 
0 60 X, 1 “X, = 0 0 0 0 
L = (J.4) 
0 0 0 0 Xz 1 | -I 0 0 
0 O 0 0 0 0 Xy 1 “Xe sal 


where the schedule terminates in flat segments represented 


by 


The optimization problem is solved by adjoining the least 
squares problem with Lagrange multipliers for the constraints 
forming the algebraic problem: 


AA OL Q Avy (3.5) 


The solution to (J.5) is determined by decomposition of the 
linear system, 

The program will optimally minimize the residuals by 
moving the breakpoints in the table. This optimization pro- 
cedure utilizes a Levenberg-Marquart nonlinear least squares 
algorithm which does not require explicit calculation of the 
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derivative of the mean residuals in each region with respect 
to the breakpoints. Also, additional breakpoints may be added 
to improve the fit. In this case, a breakpoint is inserted 

in the center of a region which has the largest mean residual 
and the optimization is repeated. 


The optimal spline fits are shown in Figure J.2 for these 
breakpoints and the mean residuals are tabulated in each region 
as shown in Table J.2. The effects of optimized breakpoints 
and the improvement in fit using additional breakpoints were 
not explored in this example. 
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APPENDIX K 


DESIGN OF FTIT ESTIMATOR 


The FTIT estimator can be designed directly from Figure 
K.1 showing the transfer functions from the two inputs to the 
filter outputs. The design procedure will assume that the 
feedback, a, is small and the gains will be calculated for 


a=0. Then, a will be chosen to provide faster step response. 


The first step is to determine the filter breaks, Ay 


ros and Az The first break, Ay> determines the frequency 


o (sec”!) 

ROOT LOCUS 

ee 
1/g(p) 


FUEL FLOW 
INPUT 


SENSOR 
SIGNAL 
INPUT 





FREQUENCY 
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Figure K.1 Asymptotes of Frequency Response of FTIT Sensor 
Compensation Showing Pole-Zero Locations as a Function 
of Estimator Gains and Feedforward Parameter 











at which fuel flow and temperature sensor inputs have the 

same effect. At lower frequencies than Ay» the sensor inputs 
have more effect on estimator output than fuel flow. The 
second root, A5> determines the bandpass of fuel flow inputs. 
The largest magnitude root, Az determines the noise attenu- 


ation frequency of both inputs. 


The first order model for the FTIT sensor is dominated by 
the slow component of the response. The estimator is designed 
to track this input to set steady performance. The higher 
frequency components of the response are determined by sched- 
ule and fuel flow inputs. The sensor model parameter, a, 
is chosen to be 1/5 sec’. 

Figure K.1 shows the root locus for the filter poles. The 
gain, kKjyy> is chosen to place the bandpass lag, An, at 
O.53= 05-5 sec! or a 2 to 3 second time constant. This will 
decouple the estimator from the transient plant roots but not 
from the integral trim roots. The value of ki is chosen to 
be 0.5 to assure this situation. The gain, kz5> determines 
the attenuation of the fuel flow channel at mid-frequency and 
the overall noise attenuation, Az. This value is chosen at 
0.5 to place thiss pole at 0s sec 1 as shown. Finally, ko 
is chosen along with ki> to set the low frqeuency break at 
A,- The value of k 


1 PAN 
mate pole/zero cancellation between a) and -a. Finally, 


is chosen at +0.1 to keep the approxi- 


choosing k,,=2.0 places A, at 0.04 sec’ 4, A, at 0.33 sec! 
and r-2 at 3.0 sec}. The resulting frequency response is 
shown in the figure. A value of a=0.1 is chosen to provide 
10:1 ratio between low and mid-frequency fuel flow response 


and an acceptable feedforward value to enhance the transient 


response. 
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